Femicide in Mexico: Statistical evidence of an increasing trend

This study analyzes whether femicide in Mexico has increased more severely than other life and bodily integrity crimes (e.g., homicide, culpable homicide, injuries, malicious injuries, abortion, and other crimes that threaten life). To achieve this, the Executive Secretariat of the National Public Security System database was cleaned and the number of femicides per 100,000 inhabitants was calculated, for the period from January 2016 to March 2022 in all states of Mexico. Through descriptive statistics, non-parametric analysis of means, and hypothesis tests, we demonstrate that the states with the highest number of femicides are the Estado de Mexico (State of Mexico), Ciudad de Mexico (Mexico City), and Veracruz; moreover, the number of femicides exhibits a growing trend while the total number of life and bodily integrity crimes does not. Finally, we forecast the number of femicides for the next five months. To our knowledge, there is no other article that analyzes the growth trend of femicide compared to other crimes. Visualizing and understanding that femicide is on the rise compared with other types of crimes can help the government and legislators generate policies that are consistent with the magnitude of the problem.


Introduction
In 1979, the UN Convention on the Elimination of All Forms of Discrimination Against Women (CEDAW) was created.As of 2001, all countries in Latin America had signed and ratified the convention [1].CEDAW spawned a series of domestic reforms within Latin America to address femicide as a criminal act distinguishable from other forms of homicide.Additionally, the Belem do Para Convention, adopted in 1994, is an exclusively Latin American human rights instrument designed to draw attention to violence against women in the region and cultivate strategies to protect women and defend their rights to live free of violence.The convention has been ratified by all but 5 states in Latin America and the Caribbean, and the five that have not yet ratified have acceded to the convention [2].
Despite international and domestic laws enshrining women's right to a life free from violence and numerous campaigns to define and decrease violence against women, levels of violence against women in Latin America are the highest in the world [3].The Mexican case is of particular interest due to the culture of "machismo" which aggravates violence against women and femicide by tacitly supporting male dominance in the home and the use of violence against partners and spouses [4][5][6].
In September 2020, Mexican feminists drew international attention to the problem of gender-based violence and femicides (the terms femicide and feminicide are used interchangeably throughout the text) in Mexico when they took over a federal human rights office [7].Mexico is recognized as having one of the highest rates of femicides internationally, a trend many suggest is increasing.Ten women per day die because of femicide in this Latin American country [8].Within academia, gender-based violence in Mexico has drawn interdisciplinary attention, attracting researchers from across the social sciences as well as urban planning engineers, public health experts, and policymakers.The presence and severity of the increase in femicides in Mexico remains a question of debate for politicians, policymakers, and the society.There is a debate in Mexico involving the current president Andres Manuel Lopez Obrador who has argued that both femicide and homicide stem from the neoliberal period, during which family disintegration and the loss of values occurred, leading to an increase in overall violence [9].Due to the importance of the problem and the controversy surrounding it, scholars are also interested in determining if an increase in Mexican femicides exists.
There are several reasons to suspect an increase in femicides in Mexico.Worldwide, femicide is increasing, a problem that was aggravated by the coronavirus pandemic.Latin America has shown particularly sharp increases in femicide.In 2020, Mexico's rate of femicide did not increase according to data from the Gender Equality Observatory, yet a shockingly high number of women were victims of the crime, with 948 out of every 100,000 Mexican women dying of femicide.Of countries in Latin America, only Brazil exhibited a higher rate of femicide (Fig 1).Furthermore, the National System of Statistical and Geographical Information (INEGI) of Mexico recently released a report demonstrating an increase in violence against women more broadly in Mexico [10].The women's movement in Mexico has also been very vocal regarding what they perceive to be an increasing trend in femicides, and the international community has voiced concern regarding the level of femicide in Mexico [11].What remains to be seen is if the data agrees with public and international impressions and overall trends in data regarding violence against women in the country, i.e., the purpose of this paper.As of 2015, 93% of defendants in domestic violence cases in Mexico were male and the majority of victims were female [13].Additionally, the percentage of women murdered as a proportion of total people murdered has been growing for more than a decade.According to Charles-Leija et al. [14], drawing on data from the Mexican government, the proportion of women murdered increased from 1.4 to 4.6 of all homicides between 2007 and 2012.This figure represents a staggering growth rate of more than 200 percent.

Literature review
Beginning in 2015, several studies discuss trends in femicides in different countries.Some, such as Frias [15], only use descriptive statistics.Karakasi et al. [16] analyze the trends in femicide in Greece for the period 2010-2021 finding that in 2021, while homicides decreased to 89 incidents per year, domestic homicides skyrocketed to 34 cases per year, reaching the highest annual number ever recorded nationally.Frias [15] argues that femicides in Mexico are increasing, with indigenous groups showing the greatest growth.Molinatti and Acosta [17] analyzed the growth trend of femicides in women who were abused in several Latin American countries, finding that in Mexico the number of femicides among women aged 15 to 59 increased between 2001 and 2011.
In more quantitatively rigorous studies, Ramos de Souza et al. [6] and Cardoso Meira et al. [18] analyze the homicide mortality of Brazilian women using negative binomial regression, the former in the last 35 years, and the latter only in the states of the Northwest region between 1980 and 2014.Torrecilla et al. [19] proposed a statistical methodology to identify temporal interferences in the count of time series and test their methodology with a study of femicides in Spain from 2007 to 2017.The authors find that femicides decreased during the period studied.Valdez Rodrı ´guez et al. [20] who use system dynamics to simulate the behavior of femicide in Mexico conclude that it will continue to grow exponentially.
Curious as to the actual number of femicides and trends in the data on violence against Mexico's women, we set out to answer two key questions.(RQ1) Have the number of femicides in Mexico increased since 2015?(RQ2) Can we attribute the increase to an overall increase in violent crimes or is the trend in femicides independent?
The rest of the article is organized as follows.In materials and methods: the data, method, and fundamental concepts are presented.In the analysis, the results were explained considering descriptive statistics and non-parametric means comparison; also, a hypothesis proportion test is performed between the years and a forecast to describe the trend; at the end of the section the discussion is displayed.Finally, the conclusion and further research is described.

Data
This analysis uses data related to the criminal incidents in Mexico from the Executive Secretariat of the National Public Security System (SESNSP) [21].The criminal incidence refers to the presumed occurrence of crimes recorded in preliminary investigations initiated or investigation files, reported by the Attorney General's Office and the Attorney General's Office of the states.
The dataset covers a period from January 2016 until March 2022 (87 months) across all 32 Mexican states.The registered crimes are divided into seven juridical types: Property crimes; Family crimes; Sexual freedom and safety crimes; Society crimes; Life and bodily integrity crimes; Personal freedom crimes; and Other affected legal assets (of common jurisdiction).

Method
To answer the research questions, the analysis was carried out in four different stages: 1. Descriptive statistics-Analysis of the number of femicides in the states, where the mean and standard deviation were obtained, analyzing the change by state in the years 2015 to 2021.
2. Comparison of whether the average femicides and the Life and bodily integrity crimes are the same during the 2015-2021 period.
3. Hypothesis test to confirm if the proportion of femicides with respect to the total of crimes is increasing.

4.
Forecasting methods to shed light on the future trend of femicides in Mexico.

Descriptive statistics
Graphs are useful for understanding data trends, however more formal analysis requires the calculation and interpretation of the numerical measures' summary, which serve to characterize the dataset and communicate the salient features.Measures of concentration include the mean and the median.
The mean for a data set x 1 , x 2 , . ..x n , which is the most widely used measure of concentration is simply the average of the data [22] and is given by Eq 1.
The median, on the other hand, is the mean value once the observations are ordered from smallest to largest.
Variability means show how far the data is from its mean; the variance measures show how much average deviation there is from the data with respect to the mean [22], as can be seen in the following equation.
The standard deviation s is the root of the variance s 2 and it has the same units of the data.

Means comparison
The simplest analysis of variance (ANOVA) is known as one-way or one-way ranking and involves analysis of data sampled from more than two numerical populations (distributions) or of data from experiments in which more than two treatments were used.The characteristic that differentiates the treatments of the populations is called the factor under study and the different treatments or populations are known as factor levels.
One-factor ANOVA focuses on the comparison of more than two populations or treatment means [23].Let I be the number of populations or treatments being compared and μ i the population mean i or the true mean response when treatment i is applied (i = 1, 2, . .., I).
The relevant hypotheses are H o : μ 1 = μ 2 = � � � = μ I against H a : at least two of the μ i are different.If H o is true if all μ i are identical, whereas H o is rejected when at least one pair is different.
If the sample sizes are assumed to be equal, let J be the number of observations in each sample, the data set is IJ.The means of each sample will be denoted by � X 1 ; � X 2 ; . . .� X I: hat is, The dot in the second subscript means that all the values of that subscript were added together while holding the value of the other subscript fixed.Likewise, the average of all observations is called the large mean.� X:: Further, let S 2 1 ; S 2 2 ; . . .; S 2 I be the sample variances, where Each observation X ij within any sample is assumed to be independent of one another.The distributions of population I are considered normal with the same variance.That is, each If the assumption of normality and equal variances between the samples are not met, other non-parametric tests such as the Kruskal-Wallis test or the Mood test can be applied [22].In the case of the Kruskal-Wallis test, let N ¼ P I i¼1 J i be the total number of observations in the data set, and suppose that all N observations are ordered from 1 (the smallest X ij ) to N (the largest X ij ).When H o : μ 1 = μ 2 = � � � = μ I is true, all observations come from the same distribution, in which case all possible assignments of the ranks 1,2,. ..N to the I samples are equally likely, and the ranks are expected to be mixed in these samples.However, if H o is false, some samples will be made up of observations that have small ranges in the pooled sample, while others will be made up of observations with large ranges.
Formally if R ij denotes the range of the X ij between the N observations and R i. and.� R i: denote, respectively, the total and the average of the ranges of the i−th sample, so when H o is true: The Kruskal-Wallis test is a measure of the magnitude to which the � R i: deviate from their common expected value Nþ1 2 , and is given by Eq 9 The values of K at least as contradictory to H o are the k values that are equal to or exceed K.This is an upper tail test, that is, the value of P = P 0 (K � k).Under H o , each possible assignment of the ranks of the I samples are equally likely.So, in theory, all assignments can be enumerated, the value of K for each determined, and the null distribution obtained by counting the number of times that each of the K values are present.When H o , is true and I = 3, J i � 6, i = 1,2,3 or I > 3, J i � 5, i = 1,2,. .., I, then K has a chi-square distribution with I − 1 degrees of freedom.This implies that a P value is approximately equal to the area under the curve w 2 IÀ 1 to the right of k.This implies that if this value of k is less than the level of significance α, H o will be rejected and at least one pair of means will be different, where reliability = 1 − α.

Two proportions hypothesis test
An individual is considered successful S in a population if it possesses some characteristic of interest.Then the proportions in two different samples can be defined as: • p 1 = proportion of successes S in a population 1.
• p 2 = proportion of successes S in a population 2.
That is, p 1 is the probability that an individual is successful in the population.Suppose that a sample size n 1 is selected in the first population and n 2 in the second.Let X be the number of successes S in the first sample and Y be the number of successes in the second.The independence of X and Y is assumed, whenever the two sample sizes are much smaller than the corresponding population, the distributions of X and Y are binomial.The natural estimator of the difference of the proportions p 1 − p 2 is the corresponding difference between the sample proportions X n 1 À Y n 2 .Devore [23] shows that if p1 ¼ X=n 1 and p2 ¼ Y=n 2 where X ~Bin(n 1 ¸p1 ) and Y ~Bin(n 2 p2 ) with X and Y as independent variables, then the expected value is so that p1 À p2 is an unbiased estimator of p 1 − p 2 and the variance is So, if the distributions of p1 ; p2 are approximately normal, the distribution of the estimator p1 À p2 can be assumed to be approximately normal.By standardizing p1 À p2 a variable Z, whose distribution is approximately standard normal, is obtained.
ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi The most general null hypothesis would be the form H o : p 1 − p 2 = Δ 0 , Δ 0 = 0 and Δ 0 6 ¼ 0 must be considered separately.The real problems are mostly of the case Δ 0 = 0, that is H o : p 1 = p 2 .When H o : p 1 = p 2 is true, let p be the combined ratio of p 1 and p 2 and = 1 − p. So, the standardized variable is: ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi This variable has approximately a normal distribution when H o is true.By replacing p and q with appropriate estimators a test statistic is obtained.If p = p 1 = p 2 instead of separate samples of size n 1 and n 2 from two different population with binomial distributions, we have a single sample of size n 1 + n 2 from a population with proportion p.The total number of individuals in this pooled sample that have the characteristic of interest X + Y.The natural estimator of p is therefore: Then p is a weighted average of the estimators p 1 and p 2 with the following formulation: Null hypothesis: H o : p 1 − p 2 = 0 Statistical value of test Z¼p 1 À p2 ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi The security test is that p1 n 1 ; q1 n 1 ; p2 n 2 ; q2 n 2 is at least 10.

Forecasting methods
There are different prediction methods; those that are used in a time series from historical values are known as extrapolation methods.In a forecasting method by extrapolation, it is assumed that past patterns and trends will continue into future months, without considering what caused the past data, they simply assume the trends and patterns will continue.
Although there are some characteristics such as trend or seasonality that help to prefer one forecast method over another, the fact is that it is accuracy measures such as MAD (Mean Absolute Deviation), MAPE and MSD that help to choose one method over another.
MAD is the measurement of the error size in units, the formula is: Y i is the real value of the series at the period i Ŷ i is the forecasting in the period i n is the number of periods MSE maximizes the error by raising the square of the differences, punishing those whose difference was higher compared to others, being suitable for periods with small deviations.
MAPE measures the deviation in percentage terms and not in units like the previous measures.It is the average of the absolute error or difference between the actual demand and the forecast, expressed as a percentage of the actual value.
Single exponential smoothing is one of the most widely used methods of forecasting because it requires little computation.This method is recommended when the data has no trend, that is, there is no cyclical variation or pronounced trend.In contrast, double exponential smoothing is a method used when the data has a trend [24].In this case two parameters are required: level and trend.The level serves to smooth the estimation of the data and the trend calculates an average growth at the end of the period [25].The formula is:

Analysis and discussion
In this section, we present the data analysis and discuss the results.

Descriptive statistics
As mentioned previously, SESNSP [21] dataset covers a period from January 2016 until March 2022 (87 months) across 32 Mexican states.The registered crimes are divided into seven juridical types: Property crimes; Family crimes; Sexual freedom and safety crimes; Society crimes; Life and bodily integrity crimes; Personal freedom crimes; and Other affected legal assets (of common jurisdiction).
During the analyzed period, 13,814,735 crimes were registered; out of these 1,809,735 cases (13.10%) belong to the Life and bodily integrity crimes (homicide, culpable homicide, injuries, malicious injuries, abortion, and other crimes that threaten life).In this juridical type of crime, 5,759 femicide cases were registered (with standard deviation 2.66), which represents 66.20 registered cases per month (per 100,000 habitants) in the whole of Mexico.The data is in the repository 10.6084/m9.figshare.22111211.Table 1 summarizes the descriptive statistics of the registered crimes for all 32 Mexican states.
Co ´digo Penal Federal [26] article 325 defines femicide as a crime when someone deprives a woman of her life for reasons of gender.Gender reasons are considered to exist when any of the following circumstances occur I.The victim shows signs of sexual violence of any kind; II.Infamous or degrading injuries or mutilations have been inflicted on the victim, before or after the deprivation of life or acts of necrophilia; III.There is a history or data of any type of violence in the family, work, or school environment, of the active subject against the victim; IV.There has been a sentimental, affective or trust relationship between the perpetrator and the victim; V.There is evidence that establishes that there were threats related to the criminal act, harassment, or injuries from the subject against the victim; VI.The victim has been held incommunicado, regardless of the time prior to the deprivation of life; VII.The body of the victim is exposed or exhibited in a public place.In the case that femicide is not proven, rules of homicide are applied.

Means comparison
We analyzed whether there was a significant difference between femicides per state during the period of study; thus, 32 samples were considered.This analysis was not possible with ANOVA because the samples did not meet the normality assumptions and the variances were  2 and 3, respectively-During the evaluated period, an overall median 1.0 case of femicides per 100,000 habitants was reported monthly.The results indicate statistically significant differences between the state medians (p < 0.000).The highest number of femicides is reported in Estado de Me ´xico with a median of 7.5 registered cases, followed by Veracruz (6), Ciudad de Me ´xico (4), Jalisco (4), and Nuevo Leo ´n (4).Even though, Ciudad de Me ´xico, Jalisco, and Nuevo Leo ´n have a median of 4 femicides, Ciudad de Mexico has 77 periods above the median, Nuevo Leo ´n 55, and Jalisco 37 (Table 2).The lowest femicide rates are observed in Baja California Sur (0, Aguascalientes (0), Campeche (0), Nayarit (0), Tamaulipas (0), Tlaxcala (0)), and Yucata ´n (0).In the second part of the analysis, the means of femicides from January 2015 to December 2021 were compared to investigate the tendency of the reported cases, the year 2022 was excluded from the analysis as the data includes only 3 months of the year.To better understand trends in the baseline levels of violence, we also examined the tendency regarding the number  Tables 4 and 5 present the results of the Kruskal-Wallis and Mood media tests for the reported cases of femicide with respect to each year.In both cases, H o is rejected (p < 0.000), i.e., there are differences between the years of the analysis.Both tests confirm the increasing tendency of femicide in Mexico, and the growth of the reported cases is statistically significant.
We analyzed whether a similar trend can be observed in the case of the reported numbers of life and bodily integrity crimes.In Tables 6 and 7, the same tests were applied.In this case, in the Kruskal-Wallis and Mood median tests, the p-value is greater than 0.05 (p = 0.097 and p = 0.070 respectively).Therefore, H o is rejected, that is, there are no statistically significant differences between the years in life and bodily integrity crimes.In other words, while in Tables 4  and 5 there are statistically significant changes in the number of femicides, in the total number of life and bodily integrity crimes no statistically significant trend is observed.Chi-Square = 87.77DF = 6 P < 0.000

Hypothesis proportion test
The tests in the previous section indicated a significant difference in the reported femicides from 2015 to 2021, and no significant difference in the case of total number of life and bodily integrity crimes during the same period.However, we felt it was necessary to assess whether femicides were growing or decreasing related to life and bodily integrity crimes in the time horizon.For this, the following hypothesis was proposed using the test statistic of Eq 13, this equation is based on the assumption of normality in the Central Limit Theorem (CLT).
According to the CLT, when n is large the probability P(a < X < b) can be calculated assuming that the distribution is normal and standardize it.
p 1 = proportion of femicides related to life and bodily integrity crimes in a year.p 2 = proportion of femicides related to life and bodily integrity crimes in the following year.Table 8 summarizes the results of the analysis.The proportion of the reported cases of femicides with respect to the total number of life and bodily integrity crimes has generally increased over the evaluated period, although it remained constant in a few periods.The growth of the registered cases of femicides was statistically higher in 2016 compared to 2015 (p < 0.000), in 2017 compared to 2016 (p = 0.005), in 2018 compared to 2017 (p < 0.000) and in 2020 compared to 2019 (p = 0.007), whereas the period from 2018-2019 and from 2020-2021 do not demonstrate statistically significant increases (p = 0.492 and p = 0.901 respectively).The tendency is clearly visible in Fig 4.
The year 2020 in Fig 4 is interesting.There are studies that discuss the 2020 situation related to the pandemic in terms of people being less vulnerable to crimes because of reduced mobility in public spaces.However, violence against women worsened considerably in this period as many women spent prolonged confinement periods with their aggressors, not only in Mexico but in other Latin American countries as well [27].

Forecasting
Finally, to make a forecast of the number of femicides that could be expected in Mexico in 2023, double exponential smoothing technique was used to analyze the time series.This approach is used when the series exhibits a trend, but not a seasonal pattern (as happens with femicides).Further, considering all the data but the decreasing weights in older observations, each month in each year is considered as a period of the time series (87 periods).Table 9 presents the forecasted expected femicide cases and the double exponential smoothing, and  It can be observed that a growth of femicide cases is to be expected.Comparing the forecast against the actual data, the forecast is even lower, but the forecast is within the confidence interval [28].

Discussion
Between 10 and 11 women are murdered every day in Mexico, while women, girls, and adolescents are kidnapped, and many are sexually violated throughout the country [29].Within the framework of International Women's Day, the National Citizen Observatory of Femicide (OCNF), a national network made up of 43 organizations and located in 23 Mexican states, warned of its concern about the prevalence of femicide in Mexico.We asked if femicide is part of a general trend of increasing violence in the country or if femicide is increasing independent of other crimes against life.Although newspaper articles show that femicide is growing [30][31][32], to our knowledge, no post-2010 research article shows that femicides are growing faster than other non-typical crimes against life as such.Additionally, we were only able to identify one scientific article that used quantitative data and techniques to describe femicide in Mexico [20].Thus, the objective of this article was to answer whether the number of femicides per 100,000 inhabitants is increasing compared to all crimes against life in Mexico.
Our results demonstrate that there is an increase in the number of femicides across most years of our analysis, except for 2018-2019 and 2020-2021.Furthermore, by employing non-  parametric means tests we found that there was a difference between the years for femicide (RQ1), but not for crimes against life (RQ2).
To assess whether femicides grew at a higher rate than the total number of crimes against life during our period of study, we carried out an analysis of proportions of the number of femicides with respect to the total number of crimes against life.We found that the proportion is decreasing in some periods and constant in others, which shows that femicides are growing at a different rate than the total number of crimes against life.In other words, the tendencies observed in femicide are categorically different than for other crimes against life and bodily harm.
The findings of our research set the foundation for researchers to examine what social, political, and economic factors are most relevant to femicide in Mexico, and to reevaluate public policies in Mexico regarding femicide.

Conclusion
Since the 90's, laws were to protect women from violence and numerous campaigns were launched to promote respect for women.Nevertheless, femicide in Mexico has increased.Mexico has one of the highest femicide rates in Latin America and worldwide.
Mexican policymakers debate whether the trend of femicide is greater than other crimes against life if femicide warrants attention as a specific and separate area of concern that policymakers should seek to directly address.To shed light on this question, we addressed, two research questions and collected data from the Executive Secretariat of the National Public Security System.This database was cleaned to describe the number of femicides and other crimes against life per 100,000 inhabitants in all states of Mexico from January 2016 to March 2022.
Using non-parametric tests (Kruskal Wallis and Mood) in addition to hypothesis tests, we found that the Mexican states with the highest number of femicides are the Estado de Me ´xico, Ciudad de Me ´xico, and Veracruz.Additionally, we found that the average number of femicides has increased since 2015 (RQ1) while that of other crimes against life has not (RQ2).
In publishing this research, we hope to draw policymakers' attention to femicide as a priority and to raise social awareness of the gravity of the problem.Furthermore, we recognize that this problem is multifactorial and plan future studies addressing the factors which impact the growth of femicides in Mexico from a multidisciplinary perspective.

0
Area under the standard normal curve to the right of the z bÞ H i : p 1 À p 2 < 0 Area under the standard normal curve to the right of the z cÞ H i : p 1 À p 2 6 ¼ 0 2 * ðarea under the standard normal curve to the right of zÞ Fig 2 presents the level of femicide across the country.
Fig 3 shows the comparison of the mean and median of both variables.The results reveal that the mean number of registered femicides has an increasing tendency, while the total crimes against life remain constant from 2015 to 2021.Femicide increased from an average of 83.53 per 100,000 inhabitants in 2015 to 181.03 cases in 2021, whereas 645 cases of life and bodily integrity crimes were reported in 2015 and 675.2 in 2021.The normality for each year of femicides and life and bodily integrity crimes was analyzed.As S2 Fig indicates, normality was not fulfilled, and the assumption of equality of variances was not met.For this reason, non-parametric tests were used to formalize the statistical analysis.

Fig 5 .
Fig 5. Forecast for femicides in the next five months.Authors' elaboration.https://doi.org/10.1371/journal.pone.0290165.g005 Is the smoothed value Y t in the time t S t 0 Is the double smoothed value Y t in the time t a t Calculates the difference between exponentially smoothed values b t Is the adjustment factor Ŷ tþm Is the forecast m periods ahead.

Table 1 . Descriptive statistics of monthly registered Life and bodily integrity crimes in Mexico per state and year.
The Kolmogorov-Smirnov test was performed to test normality, and we obtained a p < 0.001 for every sample.Thus, we concluded that the samples do not meet the assumption of normality as S1 Fig indicates.As a result, the test of variances was not performed.Two non-parametric tests were performed: the Kruskall Wallis and Mood Tests, which are shown in Tables a Incomplete data for 8 months in Oaxaca b Incomplete data for1 month in Oaxaca c only 3 months are included for all states https://doi.org/10.1371/journal.pone.0290165.t001not the same.

Table 3 . Mood median test: Femicide versus state.
of life and bodily integrity crimes.

Table 8 . Test of proportions of the number of femicides related to life and bodily integrity crimes, comparing the time horizon.
https://doi.org/10.1371/journal.pone.0290165.t008